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Abstract 

The polarisation of r's produced in Z decay is measured using 160 pb" 1 of 
data accumulated at LEP by the ALEPH detector between 1990 and 1995. The 
variation of the polarisation with polar angle yields the two parameters A e = 
0.1504±0.0068 and A T = 0.1451±0.0059 which are consistent with the hypothe- 
sis of e-r universality. Assuming universality, the value A e - T = 0.1474 ±0.0045 is 
obtained from which the effective weak mixing angle sin 2 9^ = 0.23147± 0.00057 
is derived. 
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1 Introduction 



Owing to parity violation in Z production and decay, the r leptons produced in 
the e + e~ — > Z — > t + t~ reaction are polarised. Except for very small 0(m 2 /m 2 z ) 
corrections, the r + and r~ have opposite helicities, hence opposite polarisations. In 
the present paper, P T is defined as the r~ polarisation. 

The P T dependence on the angle 9 formed by the r~ direction and the e~ beam is 
given, in the improved Born approximation at the Z peak, by 

_ , m Aril + cos 2 #) + A (2 COS 0) 
P T icosv) = 



(l + cos 2 fl) + f„4 /b (2cos#) ' v ; 

where Afb is the forward-backward charge asymmetry of r production and 

A = 2g l v g l A /[(g l v ) 2 + (g l A ) 2 ). (2) 

In the standard model, the vector (g l v ) and axial vector (g l A ) couplings of the Z to 
lepton I are independent of the lepton flavour and related to the effective weak mixing 
angle by 

^M = l-4sin 2 ^. (3) 

Therefore, the measurement of P T (cos8) allows a test of the e-r universality prediction 
At = A e and, assuming universality, gives a determination of the weak mixing angle. 

The analyses presented here use the complete set of data (160 pb -1 ) accumulated 
by the ALEPH detector at LEP I from 1990 to 1995. They supersede the previously 
published results [1-3] 

The t polarisation is extracted from five decay channels which amount to a total 
branching ratio of about 90%, namely r — > ei>T>, r — > [ii/u, r — > nu, r — > pu, r — > aiz/. 
The channels r — > iCz/ and r — > Kir°v are included in r — ► ttv and r — > pv respectively, 
and the a\ decay is observed as both a\ — > 7i + 7i^7i ± (3h) and a\ — > 7r°7r°7r ± (/i2tt°). 

The methods employed for the P T measurement in each channel are described in 
the following section. At each energy, a fit of Eq. [1] allows to determine uncorrected 
parameters A T and A e . Final corrections are needed to get the effective couplings at the 
Z, whose combinations 2g l v g l A /[(g l v ) 2 + (g l A ) 2 ] are written Ai to distinguish them from 
the uncorrected A[. They are computed with the help of the ZFITTER program 



2 Methods for r polarisation measurement 
2.1 Information from a polarised r decay 

The general method to measure the r polarisation in an optimal way || takes advantage 
of the linear dependence on the polarisation of the distribution of decay products. 

In each decay channel, the decay products are described by a set of n observables 
£ and their distribution reads 

W{t)=f(t)+P T g(£), (4) 
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with the normalization and positivity conditions ff{£)d n £ = 1, fg{$)d n £ = 0, / > 
and \g\ < f. 

With the optimal variable for the channel defined by 

the distribution in Eq. [| can be cast, with no loss of information, in the general reduced 
form 

W(u) = />) [1 + P T u] = l - [(1 + P T )W + (u) + (1 - P T )W~{u)] , (6) 

where W + and W~ are the distributions for positive and negative helicity respectively. 

For events in which both r's decay into hadrons, it is possible to get information on 
the line of flight of the two r's ||. It is shown in Ref. || that, in an ideal measurement, 
using the knowledge of the r direction to construct the u variable gives for p and a\ 
decays the same sensitivity as that for n decay (Table [l]). 

The polarisation is obtained by means of a fit of W + (u) and W~{uo) functions to 
the experimental W(u) distribution. The functions are constructed from Monte 
Carlo events and normalized to the acceptances. Non-tau background is taken into 
account in the fit. 

The validity of the standard model for the description of r decays has been checked 
elsewhere |7j by the measurement of the r decay parameters. 

In the standard model, which is assumed here, the decay distributions for the 
t — ► li/u, r — ► 7tv, and r —>■ pv channels are completely determined by Lorentz 
invariance and there is no hadronic decay model dependence in the definition of u. 
Furthermore, because the r polarisation always appears through the product %P, where 
X is the handedness of the neutrino 0, the definition of the uj variable is the same for 
r + and r~, as are its distributions when expressed in terms of the r _ polarisation 
P T = P T - = —P T +. 

The situation is more intricate in the case of the a\ decay. There, in addition to 
the decay angles, the two-pion masses are required to describe the hadronic system. 
Their dependence is embodied in four "structure functions", Wa, Wc, Wn, and We @, 
whose computation requires a model of the a\ decay. The model || that is used in 
the present study assumes a sequential a\ — > pn decay. The sign of the We function 



Table 1: Ideal sensitivities for the polarisation measurement in the r decay channels without 
and with the r direction r ||. The sensitivity is defined as S = l/Ay/W, where A is the 
statistical error expected for a sample of N events. 



Channel 


Sensitivity 




without t with f 


7TV 


0.58 0.58 


pv 


0.49 0.58 




0.45 0.58 


Ivv 


0.22 
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depends on the charge of the r so that the expression of the uj variable is not exactly 
the same for r + and r~. 

For the r — > Ivv and r — > 7rz/ decay modes, the optimal variable is the ratio 
x = -E/-Ebeam of the energy of the charged particle (Z or 7r) to the beam energy 

For the other channels, the r decay is described by one angle and the mass of the 
hadronic system. The subsequent decay of the hadron is described by two angles in 
the case of the p, and by two effective masses and three Euler angles in the case of the 
a\- If the r direction is not available, one of the angles used to describe the hadron 
decay is not measurable and the analytical form of the decay distributions is made 
more complex by the Wigner rotation between the laboratory and r rest frame helicity 
axes. The expressions used for the construction of the uj variables can be found in 
Refs. § and §. 

Although all the information relevant to the polarisation measurement is contained 
in the sole uj variable, the distributions of the other parameters are useful when 
checking the quality of the simulations and the understanding of energy calibrations. 
For example, the two angles used to describe the r — > pv decay when the r 
direction is not available are |2| cos-^r oc (2(E n ± + £ , 7r o)/£ , b ea m — 1) and cosipp oc 
(E n ± — E n o)/(E n ± + E n o). The comparison of their distributions in data and Monte 
Carlo is the best test of the quality of the measurement. 

2.2 The analyses 

Two complementary analyses have been performed, which emphasize differently the 
sensitivity of the estimators and the reduction of the systematic uncertainties. 

In the first one, the information from the opposite hemisphere of the event is used 
for background rejection when studying a r decay. This approach, which was followed 
in previous analyses [1-3], is called hereafter the "single r method". In this method, 
the information exploited to construct the uj variable comes from only one of the two 
hemispheres defined by the plane perpendicular to the thrust axis. 

As mentioned in the previous section, for events in which both r's decay into 
hadrons it is possible to get information on the line of flight of the two r's and construct 
an uj variable giving an improved sensitivity for the r —> pv and r — > a\V decays. This 
more global approach is followed in the second analysis and is called hereafter the "r 
direction method". 

These two methods use the same set of data and many common standard ALEPH 
tools. The two analyses have nevertheless been performed independently and thus are 
described in detail in Sections [| and The combination of their results is presented 
in Section §. 

3 Detector and data sets 
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3.1 The ALEPH detector 



The ALEPH detector and its performance are described elsewhere || |10| . The study of 
r polarisation requires a very special care in the alignment and calibration procedures 
to master the differences between real and simulated data at the required level. 

The charged track measurement rests on three elements, the Vertex Detector 
(VDET), the Inner Tracking Chamber (ITC), and the Time Projection Chamber 
(TPC). 

The VDET, a double-sided silicon strip detector, provides precise measurements 
close to the interaction point and, therefore, plays an important role in determining 
the r direction. 

The ITC is a drift chamber with wires parallel to the beam and a short drift 
time. Besides providing a trigger, it is an efficient tool against cosmic rays. Because 
its complete drift time is about 200 ns, to be compared to the 11 {is interval between 
bunch collisions, its sensitive time window is small. Therefore, there are few if any ITC 
hits in most triggers on cosmic rays. If present, they are not reconstructed with the 
correct timing and the reconstructed track is distorted leading to erroneous distances 
of approach to the interaction point. 

The TPC is a very large volume drift chamber. Together with the VDET and ITC, 
it brings a resolution on the track momentum which reaches 6 x 10~ 4 p (p in GeV), but 
to get this accuracy a precise calibration is needed. For these analyses, the distortions 
from the TPC are derived from the study of nonradiative events in and cos 8 

bins, where <fi and 9 are the azimuthal and polar angles of the \x~ with respect to 
the e~ beam. The departures from the corresponding observations on Monte Carlo 
are corrected for, considering that the distortions on the sum of positive and negative 
track momenta are due to field effects and that the distortions on the difference are 
due to sagitta measurement effects, i.e. alignment effects. 

Apart from momentum and angle measurements, the tracking system is used to 
measure the impact parameter relative to the beam axis, denoted do- 

The tracking detectors are surrounded by the electromagnetic calorimeter (ECAL). 
All these elements are installed inside a superconducting coil which provides a 1.5 T 
axial field. The return yoke of the magnetic field is instrumented to form a hadron 
calorimeter (HCAL). 

The ECAL is a lead/proportional wire chamber sampling device of 22 radiation 
length thickness. The insensitive regions between its modules (cracks) represent 2% of 
the barrel and 6% of the endcap areas. The anode wire signal is read plane by plane. 
The cathodes are divided in pads, making 74,000 towers pointing to the interaction 
point. These towers are read out in three sections in depth ("storeys" or "stacks"). 
The use of the energies measured with the pad towers is mandatory whenever a precise 
direction is needed but, in the case of r events, the energy collected on the wires may 
be advantageous. It does not suffer from disconnected channels ensuring that the total 
energy left in each module is properly known. This comes from the fact that the lack 
of a plane is, to first order, corrected for by the energy calibration. It has a very low 
level of noise (about 10 MeV for a complete module). The energy resolution therefore 
does not suffer from threshold effects and is slightly better than the resolution of the 
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pad signal. Another interesting feature is that the 45 plane signals provide a complete 
depth profile of the energy deposit. 

The calibration of the ECAL was performed independently in the two analyses 
presented below but the salient features are common and the technique is the same 
for pad and wire signals. First a correction map is applied which takes care of local 
gain variations related to mechanical distortions, then gain dependencies as a function 
of and 9 are applied to each module, finally the modules are inter-calibrated and 
normalized to the response observed in Monte Carlo. A problem in this procedure 
comes from the nonlinearity with energy of the calorimeter response which reaches 
4% at 50 GeV. Bhabha events provide a calibration at the beam energy but not 
a measurement of the nonlinearity which can be parametrized in first order by a 
constant a in the expression giving the true energy as a function of the measured 
one, E t = E m (l + aE m ). To extract a, clean electrons are selected and the E/p ratio 
is compared, as a function of energy, between real data and Monte Carlo where the 
saturation effects have not been simulated. This is done independently for barrel and 
endcaps. 

The HCAL is made of 23 layers of streamer tubes inserted between 5 cm thick iron 
plates. The individual tube hits are recorded digitally and the energy is read out in 
4,788 projective towers. This calorimeter is used primarily for discriminating between 
pions and muons. This does not require a specific energy calibration. For the current 
purpose, use is made essentially of the digital pattern provided by the firing of the 
streamer tubes. In particular, the pattern recognition described in Ref. is used to 
get variables describing the shape of the shower. 

The HCAL is surrounded by two double layers of streamer tubes, the muon 
chambers, providing additional information for muon identification. 

The luminosity calorimeters are not utilized in the analyses presented here. 



3.2 The data and Monte Carlo samples 



The data used in these analyses were accumulated with the ALEPH detector during 
the years 1990 to 1995 at centre-of-mass energies close to the Z mass. This set, often 
referred to as the LEP I data set, corresponds to about 160 pb _1 of data. A detailed 
description of the characteristics of each period can be found in Ref. [11]. For the 



present study no requirement on the luminosity measurement quality is made. 

Monte Carlo simulation samples were generated for the r + r~ events as well as for 
the different backgrounds through the complete detector simulation chain. A sample 
of 1.9 x 10 6 r + r _ events was generated using KORALZ For each year and energy 
of data taking a sample of about six times the real data sample was simulated with 
the same energy and the effective detector geometry of the year. 

A sample of about 4xl0 5 events was simulated using also KORALZ for 

different years and energies. About 5xl0 5 large angle Bhabha events were produced 

HI 



using UNIBAB p 



for different energies. For the two- 



and some with BABAMC 
photon interaction processes, the simulation used the generator PHOT02 |15| . Finally 
the four-fermion processes £ + £~V, defined []11| as a pair of leptons associated with a 



two particle system (V) of lower mass, were generated using FERMISV Jl6| . 
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4 Single r analysis 



4.1 Introduction 

The aim of this analysis is to reduce the systematic uncertainties, even at the price of 
some inefficiency Data have been processed separately for each year of LEP operation 
in order to take as much possible of time-variations in the behaviour of the 

detector. 

Particle and channel identification, as well as non-r background rejection, are 
achieved by means of likelihood estimators, in order to optimize the efficiencies. These 
estimators are constructed from Monte Carlo distributions of the relevant variables, 
but their efficiencies are directly measured on selected samples of data, which are used 
also to study the properties of the remaining backgrounds. 

Specific procedures have been designed for the 77 and cosmic ray background 
rejection. They are described in detail in the following. 

To cross-check the photon identification and reconstruction, an inclusive one-prong 
hadronic analysis, which infers the polarisation from the charged pion momentum 
distribution, has been devised, in addition to the standard study of the previously 
listed channels. 



4.2 Tools 



The general tools applied in the analysis are presented here. Additional information 
on their performances is given in Section |4.4| dedicated to the study of the systematic 
uncertainties. 



4.2.1 Charged particle identification 



The charged particle identification procedure (PID) described in Ref. (TT| is used for 
the non-r background rejection and for the decay channel selection. 

Charged particles are classified as electron, muon, or hadron (pions and kaons are 
treated as a single class thereafter referred to as "pions"). Since there is negligible 
probability of mistaking electron and muon (~2.5xl0 -4 ), the PID is based on two 
likelihood estimators which give the relative probability for a track to be an electron 
versus a pion (e/n), and to be a muon versus a pion {n/it), respectively. The fi/n PID 
is restricted to charged tracks with momentum above 1.3 GeV/c. An identification 
efficiency of about 99% with a rate of misidentification of about 1% is achieved for 
the /i's in the previously defined region and for the electrons inside a fiducial region 
which excludes about twice the size of the physical ECAL cracks. A study of selected 



data samples [|TT[] has shown that the efficiencies are reproduced by the Monte Carlo 
at the two per mil level. The small resultant corrections are applied in the polarisation 
measurements. 
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4.2.2 Photon reconstruction 



The reconstruction of photons from clusters of cells in the electromagnetic calorimeter 
is performed with the ALEPH standard algorithm [10] which searches for local energy 
maxima in the ECAL pad clusters. Fluctuations of a shower can generate "fake 
photons" |I7j which are artefacts of the clustering algorithm or true photons produced 
by secondary interactions in the ECAL. 

In order to reduce the number of such fake photons, an estimator of fakeness {Pf a k) 
is built from the relevant parameters |17|: the oriented distance to the closest charged 
track, the fraction of the shower energy deposited in the first ECAL stack, and the 
total ECAL shower energy. The sign of the distance is computed depending of the 
position of the shower with respect to the track bending in the r — projection. The 
reference distributions are taken from Monte Carlo, with the ECAL barrel and end-cap 
being handled independently. A cut on the fakeness likelihood estimator provides a 
partition of the photon candidate sample in two subsamples: the "fake photon" sample 
(Pfak >0.75), enriched in fake photons and the "genuine photon" sample (Pf a k <0.75), 
enriched in genuine photons. The distributions of the variables used to construct the 
Pfak estimator are shown in Fig. [I] for each of the two subsamples of photon candidates. 
A good agreement between data and r + r~ Monte Carlo can be observed. 



4.2.3 Non-r background rejection 

Non-r background in the selected samples has four main sources: Bhabha events, 
Z — > decays, Z — > qq decays, and 77 processes. As stated before, likelihood 

estimators exploiting the information of a single hemisphere are used for the suppression 
of the backgrounds that arise from the first three categories. 

The Bhabha estimator (S e+e ) provides the likelihood for one hemisphere to be 
part of a Bhabha event rather than of a t + t~. It is based on three variables: the e/ir 
probability of the track, its momentum, and the ECAL wire energy. Specific reference 
distributions are used for tracks pointing to an ECAL crack. 

The Z — > estimator ) is similar to the previous one with the j^/ir track 

probability replacing the e/n probability. 

A third estimator (S qq ) tags a Z hadronic decay in one hemisphere. It makes 
use of the opening angle between the charged tracks, the number of charged tracks, 
the number of objects reconstructed from the pattern of the energy deposited in the 
calorimeters, and the mass and energy of the jet. It is very similar to the estimator 
used for the leptonic cross section measurement and is described in Refs. |TTJ and |18j . 
Different reference distributions are taken for each ECAL region (barrel, endcaps, and 
their overlap). 

The 77 background is rejected by cuts on the three pertinent variables: visible 
energy, missing transverse momentum, and acollinearity of the two jets. The cut values 
are optimized for each final state. 

Two less important backgrounds can also contribute: £ + £~V and cosmic ray events. 
The £ + £~V events are eliminated by the same cuts that reject Bhabha, Z — > 
and Z — > qq events. 
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Figure 1: Distributions of the variables used to define estimators for genuine and fake 
photons. Data are represented by points while the histograms are the predictions of the 
simulation. Discriminant variables are described in the text. In each plot the distributions 
for data and Monte Carlo are normalized to the same number of photons. The distance d 77 
|l7fl is used in the second analysis only. 



Cosmic rays can contaminate the r — > hu and r — > [ivv channels. They appear 
as events with two collinear tracks of the same momentum coming from the vertex 
region and a low energy deposit in the electromagnetic calorimeter. A first set of 
cuts is applied to select them: ECAL energy below 2.5 GeV, \p\\ — \p^\ < 5.0GeV/c, 
min(|<iQ|, |c/q|) > 0.1cm, and \dl + d^\ < 0.5cm, where Pi and d % are respectively the 
momentum and signed impact parameter of the track i. The short gating time of the 
ITC is then used |TP[. On the sample of cosmic ray candidates selected by the preceding 



cuts, the distribution of the number of ITC hits is measured. It is found that 68% of 
the events have no hits and 88% less than eight. Since the number of ITC hits is related 
only to the arrival time of the particles, the comparison of its distributions in the two 
samples, cosmic ray candidates and r decays, is used first to estimate quantitatively 
the residual background in the r samples. Finally, the events in the previously denned 
sample are rejected, along with events satisfying the same criteria on ECAL energy 
and momenta, and N/ TC - < 8. 
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4.3 Channel selection 



The ALEPH standard leptonic pre-selection described in Ref. [IT]], which essentially 
requires a multiplicity of charged tracks between two and eight, is used first. The 
different channels are then selected in the angular acceptance (| cos#| < 0.9) with the 
help of the charged particle and photon identification. 



4.3.1 Leptonic decays 



200 




yy mass (GeV/c 2 ) Jet mass (GeV/c 2 ) 

Figure 2: (a) The 77 mass distribution for the "electron" hemispheres with a jet mass larger 
than 0.5 GeV/c 2 . The shaded area is the contribution of r — > evv taken from the t + t 
Monte Carlo, (b) The "7r ± "(7(s)) mass distribution in the muon sample before the hadronic 
veto. The shaded area shows the contribution of r —* [ivv taken from the t + t~ Monte Carlo. 
The black dots are the data, the solid (dashed) histogram the t + t~ Monte Carlo prediction 
taking (not taking) into account the correction to the hadron misidentification rate described 
m Section $~2l\ . 

The leptonic r decay sample is selected from hemispheres containing only one 
charged track coming from the interaction region and identified as electron or fi. 
In the case of electron, it is required that the charged track does not point to an 
ECAL crack. To reduce the hadronic r decay background possibly arising from particle 
misidentification, it is required that there is no genuine 7 pair with an invariant mass 
compatible with the ir° mass, and that the jet (hemisphere) mass is smaller than 0.5 
GeV/c 2 . Since bremsstrahlung is only significant for electrons, muon hemispheres with 
one or more photons are also rejected when the hemisphere mass lies between 0.5 
GeV/c 2 and 1.6 GeV/c 2 . This "hadronic veto" is illustrated by Fig. |2| which shows the 
77 mass in electron hemispheres, and the jet mass in muon hemispheres containing at 
least one photon. 

The non-r background is rejected by cuts on the above described £ e+e and £^ ^ 
estimators. To further suppress the Bhabha event contamination in the electron 
channel, the total ECAL wire energy is required to be smaller than 0.9 y/s. The 
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ECAL wire energy distribution, before this last cut, is shown in Fig. |3|. The agreement 
observed between the data and the sum of signal and background Monte Carlo shows 
the good understanding of the non-r background contamination in the electron channel. 
The two-photon background is rejected by the cuts defined in the previous section. 




ECAL- wires energy Ns 

Figure 3: Distribution of the total ECAL wire energy, normalized to the centre-of-mass 
energy, in the electron channel. The black dots represent the data, the line histogram the 
sum of all contributions, t + t~ , 77 — > e + e~, 77 — > t + t~ , and Bhabha events, given by Monte 
Carlo. Grey histograms, represent the backgrounds: dark grey for Bhabha events and light 
grey for 77 interactions. The vertical line displays the cut value. 



4.3.2 t — ► hv 

For the selection of hadronic decays without a neutral hadron, the presence of only one 
charged track coming from the vertex region, with a momentum larger than 1.3 GeV/c 
(limit of the validity of /i-PID), is required. It is also required that this track is not 
identified as an electron or a muon by the PID. Two further cuts are imposed on the 
same hemisphere: no genuine reconstructed photon and a calorimetric energy deposit 
consistent with the momentum of the charged hadron. This last cut reduces the K^n 
background. 

The rejection of the hadronic Z decays and other non-r backgrounds is performed 
by cuts on the values of the previously defined likelihood estimators computed on the 
opposite hemisphere. The two-photon background is rejected by the cuts described in 



Section 4.2.3. 



4.3.3 r — ► pv 

This selection is similar to the hp channel selection as far as the characteristics of the 
charged track and the non-r background rejection are concerned. The requirements on 
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the photons in the charged track hemisphere are opposite. The presence of one or two 
genuine photons is required. When a single photon is reconstructed, its energy must 
be greater than 0.5 GeV and a cut on the 717 mass (0.4 GeV/c 2 < < 1.6 GeV/c 2 ) is 
imposed. When two photons are reconstructed, no cut on the ti2^ mass is applied but 
the 77 mass has to be compatible with a 7r° mass (80MeV/c 2 < m 77 < 250 MeV/c 2 ). 

The quality of the simulation is illustrated in Fig. [|, which shows the 7r° (one and 
two 7(s)) energy distribution in the selected p sample along with the single photon 
fraction, and in Fig. |5|, where the jet mass in the one 7 sample (before the above 
described cut on this mass) is displayed. 
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Figure 4: (a) Distribution of the 7r° energy data (black dots) and t + t~ Monte Carlo (empty 
squares), (b) Data over Monte Carlo ratio, (c) Single photon fraction in the 7r° sample as a 
function of the 7T° energy. 



4.3.4 r — ► a x v, a x — > h2n° 

Here also, the selection criteria on the charged track and the non-r background rejection 
are the same as in the hu case. The additional requirements on the photons in the 
selected hemisphere are as follows. The number of genuine reconstructed photons, n T , 
ranges from two to four. For n 7 = 2, the 77 mass has to be larger than the it mass and 
inconsistent with it. For n 7 = 3, the reconstruction of at least one tt° is required and, 
for n 7 = 4, the reconstruction of two 7r°'s is required. The choice of the best association 



of the photons and the ir° reconstruction are performed by a kinematic fit [nj which 
also improves the energy determination. 

The purity of the selection and the quality of the reconstruction for the four photon 
sample are shown in Fig. ^| where the mass of a 77 system is displayed when the mass 
of the two remaining photons is consistent with the tt° hypothesis. 
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Figure 5: Distribution of the ?r 7(3) mass for data (black dots) and r + r~ Monte Carlo 
(histogram), in the r — ► sample. The shaded histogram shows the contribution of r 
background. 




Figure 6: Distribution of the 77 mass in four photon hemispheres when the effective mass 
of the two other photons is in the ir° region. Data are represented by black dots and Monte 
Carlo by histograms. The shaded histogram shows the contribution of r background. 
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4.3.5 r — ► a^v, — » 3h 

The presence of three charged tracks coming from the interaction region is required. 
At least two of the three tracks must be identified as charged hadrons by the PID. The 
following cuts against the presence of 7r°'s in the hemisphere are applied: no genuine 
photon with energy larger than 1 GeV, no genuine 7 pair with > 3 GeV, and 
m 77 compatible with m n o. Hemispheres containing a V°, as defined by the ALEPH 



V° package |I0| , whose mass is compatible with a K s are rejected. In addition, the 
effective mass of at least one of the 7r~7r + combinations is required to be larger than 
0.55GeV/c 2 . The non-r background rejection is the same as in the other hadronic 
decay channels. 



4.3.6 Inclusive channel r — ► h > 07^ 

The charged track selection and the non-r background rejection are those used in the 
selection of the hp channel. Only the cuts against the two-photon event background 
are tighter, because the pion momentum distribution used for the polarisation 
measurement is very sensitive to this background. Since the purpose of the study 
of the inclusive channel is to cross-check the photon and 7r° treatment, no cut using 
information from the electromagnetic calorimeter in the selected hemisphere is applied. 

4.4 Systematic uncertainties 

The sources of systematic errors affecting the r polarisation measurement are: 

• r selection efficiency, 

• non-r background contamination, 

• 7/71"° reconstruction, 

• energy /momentum calibration, 

• modelling of a.\ decay and radiative processes, 

• r decay channel cross-talk. 

A summary of the errors in the single-r method is given in Table ^| and their estimation 
is discussed in the following. 

The dominant contributions are photon/7r° reconstruction, ECAL calibration, non- 
r background, and modelling of the decay for the a%. This last effect is evaluated @ 
by changing the model and its parameters within the limits allowed by the data. 

The systematic effects are evaluated for each of the previously described tools using 
the difference between the data and Monte Carlo plus one standard deviation as input. 
The errors are propagated through all the steps of the analysis, taking into account the 
fact that a given tool can be used in several places. For that, A T and A e are measured by 
means of the standard W /± (co | ) functions on Monte Carlo events weighted for the change 
of the tool response. The shifts of the values of A T and A e give the associated errors. 
Cross-checks of the estimations are presented for some of the important contributions. 
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Table 2: Summary of the systematic uncertainties (%) on A T and A e in the single-r analysis. 



Source 


h 


P 


3h 


h2n° 


e 


P 


Incl. h 


selection 


- 


0.01 


- 


- 


0.14 


0.02 


0.08 


tracking 


0.06 


- 


0.22 


- 


- 


0.10 


- 


ECAL scale 


0.15 


0.11 


0.21 


1.10 


0.47 


- 


- 


PID 


0.15 


0.06 


0.04 


0.01 


0.07 


0.07 


0.18 


misid. 


0.05 








0.08 


0.03 


0.05 


nhotou 


0.22 


0.24 


0.37 


0.22 










n 1 q 


u.uo 


u.uo 


U.IO 


n £\A 

U.O^i 


U.U 1 


n 1 ^ 

U. 1U 


t RR 


n dq 

U .UV 


n fu 


n 1 n 

U. 1U 


u. 


u.uo 


u.uo 


U. 1 o 


modelling 

111V vivlllllk-. 






0.70 


0.70 






0.09 


MC stat 


0.30 


0.26 


0.49 


0.63 


0.61 


0.63 


0.26 


TOTAL 


0.49 


0.38 


1.00 


1.52 


0.96 


0.93 


0.87 








A e 










Source 


h 


P 


oh 


h2ir° 


e 


P 


Incl. h 


tracking 


0.04 










0.05 




non-r back. 


0.11 


0.09 


0.04 


0.22 


0.91 


0.24 


0.17 


modelling 






0.40 


0.40 








TOTAL 


0.12 


0.09 


0.40 


0.47 


0.91 


0.25 


0.17 



4.4.1 Charged particle identification 

The comparison of the PID performances on data and Monte Carlo is reported in 



Ref. ||1 1|| . The uncertainty on the data/Monte Carlo efficiency-ratio computed for each 



year of data taking is used to estimate systematic error. 



4.4.2 Photon reconstruction 



Clustering algorithm 

The first step of photon reconstruction is the "standard " ALEPH clustering 
algorithm |T{J. Its parameters are the distance to the closest charged track and the 



threshold energy. It has been shown in Ref. [17] that the distance parameter can 
be changed by 1 mm, while preserving within one standard deviation the agreement 
between data and Monte Carlo on the distance distribution. This parameter shift 
introduced in the simulation induces negligible effects on the polarisation, even in 
channels containing photon(s). The possible effect of the threshold value is investigated 
by modifying the reconstructed shower energy in the Monte Carlo by the clustering 
correction |T0]. The resulting variation of the polarisation measurement is small but 



not completely negligible. It is part of the value given in the line "photon" of Table |2|. 

Fake photon estimator 

The second step is the genuine/fake classification for which the Pf a k estimator 
is used. In the first place, the efficiency of the identification of genuine photons 
is measured on a selected sample made of photons forming a 7r° with one "high 
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quality" photon (converted 7 or isolated shower with 4 GeV< E < 20 GeV). For energies 
greater than 4 GeV, the efficiency is not significantly different from one. The efficiencies 
for E < 4 GeV are shown in Fig. [7]. The ratio between data and Monte Carlo is 
(99.50±0.19)x 10~ 2 and an efficiency variation of 0.69xl0~ 2 is used as input in the 
error computation. 
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Figure 7: (a) Efficiency of the fake photon estimator for genuine photons. The black dots 
are the data, the empty squares the t + t - Monte Carlo, (b) Data/Monte Carlo ratio for this 
efficiency. 

Since the Monte Carlo is able to reproduce the characteristics of the fake photon 
production in a hadronic shower but not its precise normalization fi~7| , a weighting 
procedure is used to correct it. The weight for an event is written w(p) n , where n 
is the fake photon multiplicity and p the hadron momentum. The function w(p) is 
obtained by a fit of the likelihood estimator distributions. A check of the procedure is 
given by the comparison of data and Monte Carlo for the 77 mass spectra shown in 
Fig. |[ The related systematic error is estimated in the standard way described at the 



beginning of Section |4.4| , comparing the measured weight to one. 



4.4.3 Bhabha, , and qq event backgrounds 

The 8 e+e , , and 8 qq estimators use information from only one of the 

hemispheres of an event, so correlations with the other hemisphere can only be due 
to geometrical correlations in the detector construction (ECAL cracks, overlap, etc.). 
The measurement of the estimator efficiencies in the non-correlated regions is described 
first, and then an evaluation of the residual background, which is not affected by the 
possible correlations, is presented. This last point concerns only the hadronic channels 
because the selection criteria for leptonic r decays eliminate the correlated, insensitive 
parts of the detector. 

Efficiency of the estimators 

In order to measure the efficiencies, almost pure samples of data and Monte Carlo 
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Figure 8: The 77 mass spectrum, for a sample of ^77 hemispheres: (a) when the two 
photons are classified as fake, (b) when 1 photon is classified as fake and the other as genuine. 
The black dots are the data, the dashed line histogram the t + t~ Monte Carlo and the full line 
histogram the t + t~ Monte Carlo, with a weighted fake photon contribution. The hatched 
histogram is the fake photon contribution in the t + t~ Monte Carlo, without weighting (see 
text). 



events are selected for the t + t~, e + e~, and qq final states. A clean r decay 

(three charged 7r's or a p), a high energy electron, a high energy muon or a hadronic 
jet is required in one hemisphere and the opposite one is used to study the estimators. 
Examples of their responses in the case of S e+e and E^ +tJ " are shown in Fig. B] and [TO. 



A cut at a value of 0.8 on these two estimators allows a good separation between 
r's and leptonic background. The measured efficiencies, taking into account the small 
residual contamination of the samples, are given in Table [| These numbers are used 
to determine the systematic uncertainties which contribute to the line "selection" of 



Table 3: Efficiencies of the £ e e and f ^ estimators (in %) measured on data and Monte 
Carlo samples 







"e 




DATA 


M.C. 




98.92 ±0.03 


98.95 ±0.01 


e + e~ 


2.77 ±0.02 


2.66 ±0.02 










DATA 


M.C. 


T + T - 


99.45 ±0.03 


99.48 ±0.01 




0.17 ±0.04 


0.16 ±0.05 
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Figure 9: (a) The S e e distribution after tagging an electron in the opposite hemisphere, 
(b) Efficiency on Bhabha events as function of the cut value. Black dots are the data, empty 
boxes and histogram, the Bhabha Monte Carlo. 
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Figure 10: (a) The f M distribution after tagging a r in the opposite hemisphere, 
(b) Efficiency on t + t~ events as function of the cut value. Black dots are the data, empty 
boxes and histogram, the r + r~ Monte Carlo. 
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Table |2| and to the line "non-r background" for the two leptonic channels. 



A similar analysis, already described in Ref. |TTJ, is performed for the S qq estimator. 

Residual backgrounds 

The estimation of the residual background, chiefly Bhabha events, in the hadronic 
channels is an essential task on account of its impact on the A e measurement. For that 
reason, procedures were developed to measure directly the background contributions 
from the data themselves. 

The method is exemplified by the case of the r — > hv channel. Applying all the cuts 
used for the channel selection but the e/ir identification cut on the r side, and the S e+e 
cut on the opposite one, a set of data is constituted which is divided in two subsets 
by the response of the e/ir PID: (a) for e PID and (b) for 7r PID. The distributions 



of £ e e in the two subsets and the r + r Monte Carlo expectations shown in Fig. [11 
clearly display the Bhabha event background. 




Figure 11: The distribution of the £ e e estimator in the hv channel, (a) for events with 
an electron identification, (b) for a pion identification. The shaded histograms are the t + t~ 
Monte Carlo predictions. 



The events with £ e+e > 0.8 in sample (b) are used to determine the kinematic 
characteristics of this background. The normalization of the residual background 
after the 8 e+e cut is given by the number of events with S e+e > 0.8 in sample (b) 
multiplied by the ratio of the numbers of events with S e+e < 0.8 and S e+e > 0.8 in 
sample (a). This background is subtracted for the polarisation measurement and the 
statistical error on its estimation enters the systematic error. The small difference of 
the efficiencies of the cut on S e+e in samples (a) and (b), mainly due to geometrical 
effects and initial state radiation, can be computed by Monte Carlo and is added to 
the systematic uncertainty. 

As a cross-check, an independent study of the Bhabha event background has been 
performed, using the tools developed for the measurement of the r leptonic branching 
ratios |T9 |. Its conclusions are consistent with the above presented results. 
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Similar procedures are used to study the contamination. They show that 

such a background is reduced to a negligible level by the selection process. 



4.4.4 Other backgrounds 

Two-photon background 

In a similar way, the residual 77 background can be studied in the data from the 
excesses of events in the distributions of the variables (acollinearity, visible energy, etc.) 
used for its rejection. 

£ + £~V background 

Such a background appears mostly in the leptonic channels. The contaminations 
evaluated by a study of the opposite hemispheres are (0.2 ± 3.2) x 10~ 4 for the 
electron channel and (5.5 ± 3.6) x 10 -4 for the muon. To compute the tiny associated 
uncertainties the FERMISV Monte Carlo fl6|] is used. 



Cosmic ray background 

The contaminations of the two relevant channels, r — > hv and r — > \wv are 
estimated using the independent information from the track parameters and the 
number of ITC hits. Their values are (6.9±1.9) x 10 -4 and (1.5±0.5) x 10~ 4 respectively, 
leading to very small contributions to the uncertainties. Out-of-time cosmic rays are 
normally identified as hadrons, owing to the HCAL and muon chamber inefficiency in 
such events ||19| . 



r background 

Uncertainties due to the crosstalk between different decay channels are evaluated 
by varying within their errors the different branching ratios used in the Monte Carlo. 



4.4.5 Detector calibration 



Momentum calibration 



The uncertainty on the momentum calibration described in Section |3J] is used, for 
each cos 9 bin, to estimate the contribution to the polarisation systematic errors. 

ECAL calibration 

The ECAL energy calibration (Section |3.1|) is obtained from electrons taken in 
Bhabha events, 77 processes, or r decays to electrons. Uncertainties coming from the 
statistics of these samples and the handling of the ECAL saturation are introduced as 
systematic errors, leading to an uncertainty of 2.5 x 1CU 3 on the calibration, which is 
used as input for the computation of polarisation errors. 

An independent check of the ECAL calibration has been performed using the /i + /i~7 
events for which the photon energy is deduced, with a small error, from the angles, by 
the formula 



E 



<ni, . _ sin0 M+M - v '~ 



sin + sin 9 ll + 1 + sin 9 ltl - 



The distributions of AE 7 /a^, where AE^ is the difference £^ alc - — E™^ between the 
kinematically calculated and ECAL energies, and cr^ the expected resolution, are 
shown in Fig. ^ for data and Monte Carlo and a very good agreement can be observed. 
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Since some calorimetric cuts are used in the selection of the two hadronic final states 
without photons, a tuning of the ECAL response to charged hadrons is also performed. 
The related uncertainties are independently evaluated. 




-10 o 10 20 30 
A(Ey)/o(Ey) 



Figure 12: Difference between the reconstructed and calculated photon energy, AE^/aEf 
for ^ + /i _ 7 events. The insert is a zoom of the region [-5,5]. The black dots represent the 
data and the empty squares the /U + //~7 Monte Carlo. 



4.5 Results 

Table [| gives the numbers of selected data events, the selection efficiencies in the 
angular acceptance, and the contaminations from r and non-r backgrounds. 

The polarisation is measured independently for each channel and for each year of 
data taking. Some of the systematic errors are common to all the years and some are 
common to different channels. 

In a first step, the measurements for each channel in the different LEP operation 
periods are combined, taking care of the common systematic errors and adding the 



small corrections dependent on the centre-of-mass energy described in Section |0 
The values obtained are given in the corresponding lines of Table 



Figure [13] shows the x distributions for the hadron, muon, and electron channels; 
and the u distributions for the pv and a\v channels (3h and h2ir°). The distributions 
for data and the contributions (Monte Carlo) from the two r helicities and non-r 
background are superimposed. 

The measurements from all the channels are then combined taking into account 
their correlations evaluated by Monte Carlo. The inclusive hadronic analysis has been 
devised as a cross-check of the exclusive ones since their systematic uncertainties 
are essentially uncorrelated. The values given in Table [5] show that this test is 
satisfactory. Statistically, on the contrary, the exclusive and inclusive analyses are 
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Figure 13: Distributions of the x variable for r — > evv (a), r — ► /ii^z/ (b), and r — > hu (c). 
Distributions of the u> variable for r — > pv (d) and r — > ai^: ai — > 3/i (e) and ai ^> /7,27r° (f). 
The dotted and dash-dotted lines correspond to the contributions of left- and right-handed 
r's respectively. The shaded area is the non-r background contribution. 
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Table 4: Number of selected events in the data, selection efficiency within the angular 
acceptance, and contamination from the r background and the non-r background. In the 3h 
mode, the KKir and Ktttt modes are treated as r background. 



Channel 


h 


P 


3h 


h2n° 


candidates 


34249 


75296 


27854 


27924 


efficiency (%) 


76.6 


78.1 


77.2 


66.9 


rback. (%) 


10.6 


10.9 


10.9 


24.3 


Bhabhaback. (%) 


0.16 ±0.03 


0.07 ±0.02 


0.06 ±0.05 


0.10 ±0.08 


other non-r back. (%) 


0.24 


0.11 


0.11 


0.09 



Channel 


e 


P 


Incl. h 


candidates 


47106 


50585 


134108 


efficiency (%) 


75.9 


86.4 


63.3 


rback. (%) 


0.8 


0.7 


0.5 


Bhabha back. (%) 


1.00 ±0.03 




0.07 ±0.01 


other non-r back. (%) 


0.50 


0.30 


0.40 



Table 5: Results for A T and A e obtained with the single-r analysis. The first error is 
statistical, the second systematic. 



Channel A T (%) A e (%) 



hadron 


15.21 


± 


0.98 


± 


0.49 


15.28 


± 


1.30 


± 


0.12 


rho 


13.79 


± 


0.84 


± 


0.38 


14.66 


± 


1.12 


± 


0.09 


al(3h) 


14.77 


± 


1.60 


± 


1.00 


13.58 


± 


2.11 


± 


0.40 


al(h27r°) 


16.34 


± 


2.06 


± 


1.52 


15.62 


± 


2.72 


± 


0.47 


electron 


13.64 


± 


2.33 


± 


0.96 


14.09 


± 


3.17 


± 


0.91 


muon 


13.64 


± 


2.09 


± 


0.93 


11.77 


± 


2.77 


± 


0.25 


pion inclusive 


14.93 


± 


0.83 


± 


0.87 


14.91 


± 


1.11 


± 


0.17 


Combined 


14.44 


± 


0.55 


± 


0.27 


14.58 


± 


0.73 


± 


0.10 
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strongly correlated so that the inclusive measurement brings only little new information 
in the final result. 

The final values of A e and A T are given in the last line of Table [|. Their errors are 
corrected for the correlation between opposite hemispheres. 



As an illustration of the understanding of the p channel, Fig. |bj shows the r and 
p decay angle distributions with good agreement between data and Monte Carlo, even 
for low energy p's (cosVv < —0.9). 



5 Analysis with the r direction method 
5.1 Philosophy of the method 

The t polarisation analysis described in this section is based on two ideas: use 
of proven analysis tools and enhancement of the measurement sensitivity. On one 
hand, event selection, particle identification, photon and 7r° reconstruction, and decay 
classification are taken from previous ALEPH analyses of r leptonic and hadronic 
branching ratios [ 17] , [131 , and of the hadronic spectral functions |2(], |21j. On the other 
hand, the maximum sensitivity to the r polarisation can be achieved in the hadronic 
decay channels through a complete set of observables including the determination of the 
initial r direction. This new information, as compared to standard analyses, is obtained 
through a kinematic reconstruction complemented by the precision measurement of the 
charged particles using the microvertex detector in order to remove part of the intrinsic 
ambiguity of the method. 



5.2 Event selection and decay classification 
5.2.1 r-pair selection 

This analysis is based on a global event selection that retains r pair candidates from 
Z decays. Each event is divided into two hemispheres along the thrust axis. The r 
event selection is described in 11171 and references therein. The method is based on the 



removal of non-r background identified as Bhabha, mu-pair, hadronic, photon-photon 
induced, and cosmic ray events. A well-tested procedure of reconstructing the spatial 
energy flow 



10 



is used in order to include all the energy measurements and their 
topology in the detector. 

However, to enhance the selection efficiency and to reduce the biases in the 
polarisation measurement, 
introduced: 



some minor modifications of the usual cuts [IT? 



arc 



• cuts against 77 events and the cut on the energy of the leading tracks in both 
hemispheres against Bhabha and di-muon events are not applied when both 
hemispheres are identified to be hadronic. 

• in order to reduce the contamination from 77 — > rr and 77 — > qq background, 
events with one hemisphere with a single charged hadron and no 7r° are removed 
if I cosfll > 0.792 and x < 0.088. 
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Figure 14: Distributions of the two decay angles in the r — > pv channel. The vertical axes 
are in unit of events/0.1. The dotted and dash-dotted lines correspond to contributions of 
left- and right-handed r's respectively. The non-r background contribution is negligible. 
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hadronic cuts are loosened taking advantage of the more precise hadronic mass 
determination after ir° reconstruction, rather than relying on the energy-flow 
variables used in the selection. 



5.2.2 Particle identification 

A likelihood method for charged particle identification incorporates the relevant 
information from the detector. In this way, each charged particle is assigned a set 
of probabilities from which a particle type is chosen. A detailed description of the 
method can be found in Refs. |19| and |f22 |. 



Briefly, the particle identification is based on the following information: (i) the 
dE/dx in the TPC, (ii) the longitudinal and transverse shower profile in the ECAL near 
the extrapolated track, and (in) the energy and the average shower width in the HCAL, 
together with the number of planes, out of the last ten of the HCAL, which had fired and 
the number of hits in the muon chambers. Because of the calorimetric requirements, 
the complete identification method is only applied to particles with momenta larger 
than 2 GeV/c. 

The performance of the identification method has been studied in detail using 
samples of leptons from Bhabha, /i-pair, and two-photon events, and hadrons from 
7r°-tagged hadronic r decays. All data sets cover the full angular and momentum 
range |17|, |19| . Ratios of particle identification efficiencies in the data and the 
corresponding Monte Carlo samples are obtained as a function of momentum. They 
are then used to correct the Monte Carlo efficiencies in the r + r~ sample, thus 
taking into account the different momentum distributions and particle environment. 
Misidentification probabilities are obtained in the same way. The results of these 
analyses were already given for the 1990-1993 data [Tj|. Figures 15 and W show the 



corresponding efficiencies and misidentification probabilities measured in the 1994-1995 
data. 

Table |6] gives the efficiency matrix for one-prong r decays used for this analysis, i.e., 
for particles with a momentum larger than 2 GeV/c and not in a crack region between 
ECAL modules. The 2 GeV/c cut is not applied to electron candidates because of 
the good dE/dx separation at low momentum, while the crack veto is not required 
for muons which are essentially identified in the HCAL. The numbers refer to the 
1991-1995 data sample. 



Table 6: Particle identification efficiencies and misidentification probabilities (in percent) in 
one-prong r decays as corrected from control data samples. The values are integrated on the 
corresponding momentum spectra (h refers to the single hadron channel) above 2 GeV/c, 
excluding the cracks between ECAL modules. 



Id. X True — > e /j, h 

99.57 ± 0.07 <0.01 0.71 ± 0.04 

<0.01 99.11 ± 0.08 0.72 ± 0.04 

0.32 ± 0.06 0.88 ± 0.08 98.45 ± 0.06 

0.09 ± 0.04 0.01 ± 0.01 0.12 ± 0.03 



e 
h 

no id. 
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Figure 15: The identification efficiencies for electrons and muons as a function of momentum 
for 1994-1995 data (in percent). The efficiencies simulated by the Monte Carlo for the r 
sample are corrected by measurements performed on control samples as discussed in the text. 
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Figure 16: The identification efficiency and misidentification probabilities for hadrons as a 
function of momentum for 1994-1995 data (in percent). The values are obtained directly 
from data using a 7r°-tagged hadron sample. 
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5.2.3 Non-r background 

A new method is developed to directly measure in the final data samples the 
contributions of the major non-r backgrounds: Bhabhas, fi + fi~ pairs, and 77 — > 
e + e~,/i + /i~ events. The procedure does not rely on an absolute normalization for 
the simulated channels, which is difficult to obtain reliably after the large necessary 
rejection is applied. It only makes use of a qualitative description of the distribution 
of the discriminating variables. The basic idea is to apply cuts on the data in order to 
reduce as much as possible the rr population while keeping a high efficiency for the 
background source under study. The background fraction is directly fitted in the cut 
data sample and then extrapolated to the standard sample. 

For Bhabha background, different event topologies are considered in turn. For e-e 
and e-crack topologies, the acoplanarity angle (defined by the two tracks azimuthal 
difference) is required to be larger than 179° or the acollinearity angle (i. e. 180° 
minus the angle between the directions of the two tracks) should be less than 175° 
with the difference of transverse energies being less than 3 GeV. The corresponding 
Bhabha efficiency (compared to the normal sample) is 88% as determined from the 
simulation. The angular distribution of the restricted sample is then fitted to rr 
and Bhabha components (also including a small contribution from the other non-r 
backgrounds) from the simulation. Therefore, the Bhabha Monte Carlo input is only 
used to determine the (large) selection efficiency and the cos# shape inside the final 
sample, not relying on any determination of the absolute Monte Carlo normalization. 
The derived Bhabha contribution has a statistical uncertainty which is assigned as a 
systematic error. Several combinations of variables have been tried, showing a good 
stability of the result within its error. For the more numerous e-h and h-h samples, a 
total energy cut of 55 GeV is applied and the main reduction of r events is achieved 
by suppressing true hadrons as compared to electrons misidentified as hadrons: this is 
done by restricting the opposite hadron in an e-h event to have an electron identification 
probability larger than 0.01 (most of the true hadrons are below this value). A similar 
approach is used for the h-h samples with even tighter identification probability cuts 
for hemispheres with only one charged hadron (in this case 86% of the hadrons have 
an electron probability smaller than 0.0002). The Bhabha signal then appears on a 
distribution of the acollinearity angle between the two hemispheres (89) at values close 
to 180°. 

A similar technique is used to estimate the /i-pair background. An efficiency of 
90% is obtained by requiring an acoplanarity angle larger than 178° for events in 
the final sample. Here, the fitted distribution is that of the calculated photon energy 
along the beam for a postulated /i/i77beam kinematics to take the most general case 
compatible with only the information on the two muons. For [i-h and h-h final states, 
the procedure follows the one used for Bhabha background. 

The measurement of the remaining background induced by 77 interactions is 
based on cuts on the total event energy (typically E < 35 GeV, depending on the 
channels) and the acollinearity angle between the two hemispheres (89 < 175°). The 
77 background shows up clearly at small overall transverse momentum well above the 
small r surviving contribution. 

All the above non-r background contributions are determined on the polarisation 
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data sample with the help of special cuts. The corresponding Monte Carlo generators 
are only used to compute the efficiency of the cuts. This method is illustrated in Fig. [17 
showing the evidence for Bhabha and 77 — ► rr or qq backgrounds in the pv channel, 
the contribution from pp being negligible in this case. 




Figure 17: Distributions in the final p sample after applying cuts to reduce the r contribution 
and enhance non-r background (see text for details) : (a) acollinearity distribution to enhance 
Bhabha background; (b) acollinearity distribution to enhance pp background; (c) total 
transverse momentum normalized to the centre-of-mass energy to enhance 77 — > hadrons 
background. The histograms give the respective r contributions as expected from the 
simulation. 



As an important bonus, this procedure leads to a determination from data alone 
of the background distributions which are relevant for the polarisation measurement. 
This method is preferable to using the Monte Carlo predictions which are questionable 
in view of the large rejection factors achieved. Distributions of the u variable obtained 
in the pv channel are given in Fig. [18|. 

As in the analysis of Section f|, cosmic ray background is determined to be very 
small after the additional cuts ||19|| . The contamination from hadronic Z decays affects 
only the hadronic channels. Its small contribution is estimated with the JETSET 7.4 
generator [24]. Tests were previously made to ascertain the reliability of the prediction 
for the rate of surviving low multiplicity events by comparing to data |T7J. A 30% 
uncertainty was derived from these tests. 
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Figure 18: Distributions of the to variable in the final p sample after applying cuts to 
reduce the r contribution and enhance non-r background (see text for details): (a) Bhabha 
background; (b) background from 77 — * hadrons. In both cases the small remaining r 
contribution estimated from the simulation has been subtracted out. 

The non-r background contributions are given in Table |7|, separately for the 
considered channels. 



5.2.4 Photon and 7r° reconstruction 

As in the previous analysis, the photon and 7r° reconstruction is performed according to 
a likelihood method in which several estimators and probabilities are computed in order 
to distinguish between genuine and fake photons produced by hadronic interactions in 
ECAL or by electromagnetic shower fluctuations Jl7| . 



An algorithm is used to form photon pairs in each hemisphere |T7| using a 



7T 



estimator Dfj 



Py.P^.P^ for two photons i and j, where P n o is the probability from 
a kinematic 7r°-mass constrained fit. A good agreement between data and Monte Carlo 
is observed for the distributions of P n o and DJj . Figure ^ shows the corresponding 
7T° mass distributions with the expected contribution from photon pairs including at 
least one fake photon. 

High-energy 7r°'s with overlapping photon showers are reconstructed through an 
analysis of the spatial energy deposition in the ECAL towers ("unresolved 7r°'s"). 
All the remaining single ("residual") photons are then considered and those with 
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Table 7: Channel identification efficiencies and contaminations as measured on the data 
for the r direction method. The efficiency includes a geometrical acceptance factor of 0.86. 
"acol" refers to the sample of events with a leptonic r decay, where information from the 
acollinearity angle || |2j| is used. 



Channel 


h 


P 


3h 


1 r\ n 


candidates 


33350 


78553 


25287 


28757 


efficiency (%) 


63.1 


69.1 


62.9 


57.5 


rback. (%) 


7.3 


9.1 


3.9 


20.7 


Bhabha 


57 ± 10 


118 ± 40 




9 ± 5 


/i/i 


20 ± 5 


16 ± 14 




2 ± 2 


qq 


61 ± 30 


205 ± 69 


72 ± 35 


37 ± 15 


jjee 








2 ± 3 




7 ± 7 








jjtt + jjqq 


15 ± 10 


31 ± 11 


± 4 


3 ± 5 



Channel 


e 


n 


acol 


candidates 


52952 


50249 


85035 


efficiency (%) 


68.7 


74.8 


74.4 


rback. (%) 


0.7 


1.0 


0.7 


Bhabha 


383 ± 66 




187 ± 37 


mi 


10 ± 3 


212 ± 55 


94 ± 24 


qq 


8 ± 4 




8 ± 4 




187 ± 45 




91 ± 20 




7 ± 5 


241 ± 57 


106 ± 24 


jjtt + jjqq 


22 ± 6 


22 ± 6 


28 ± 7 



P 7 > 0.5 are selected as 7r° candidates. The P 7 distributions are given in Fig. |2^ for all 
photons and for the residual ones after renormalizing the amount of fake photons in the 
simulation by a factor depending on the r decay hadronic final state and ranging from 
1.17 to 1.28. The shapes in data and Monte Carlo are in good agreement. The energy 
dependence of the data distributions is also well reproduced by the simulation: hence 
the photon efficiencies resulting from the D^- and P 7 cuts are adequately described by 
the Monte Carlo, as already observed in Section |4.4.2| and Fig. [7|. 

Finally, photon conversions are identified following the procedure described in 
Ref. [I7Q. They are added to the list of good photons and included in the ir° 



reconstruction. The fractions of resolved, unresolved and residual 7r°'s are presented 
in Fig. |21] as a function of the ir° energy. It is seen that the data distributions are well 
reproduced by the simulation. The fraction of resolved "calorimetric" 7r°'s remains 
significant out to large energies in this analysis as compared to the previous one (see 
Section 4.3.3 and Fig. |]); this effect is the result of a looser 7r° definition here, where 
highly asymmetric 77 pairs are more likely to be counted as resolved 7r 0, s, whereas they 
would be identified as unresolved in the other method because of the tight invariant 
mass cut. 
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Figure 19: Photon-pair mass distribution for reconstructed 7r° candidates. Data (points) 
are compared to the simulation (histogram). The shaded histogram gives the simulated 
contribution where at least one of the photons is fake. 




Figure 20: Photon estimators used to distinguish between genuine and fake photons (P 7 = 1 
corresponds to the case of genuine photons): (a) estimators for photons used in resolved 7r° 
reconstruction (the photon energy is not used), (b) estimators for residual photons (using the 
photon energy). Data (points) are compared to the simulation (histogram). 
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Figure 21: Fractions of different 7r° types as a function of the 7r° energy. The points represent 
the data and the open squares the simulation. Also indicated with stars are the fractions for 
resolved tt° containing at least one converted photon. 



5.2.5 Decay classification 

r decays are classified according to the number of charged tracks and their 
identification, and the number of reconstructed 7r°'s. Beyond the classification achieved 
in the previous analyses [17,19-21], some additional cuts have been developed in order 
to reduce the p and K* background level in the h channel, as explained below. 

In order to reject p decays where the 7r° is close to the charged pion and undetected 
as such, the sum of the energies in the first two stacks of ECAL is required to be less 
than 75% of the charged pion momentum. The energy is measured including all the 
pads in a cone of 30° around the track. To further reduce the p contamination, the 
hemispheres with one charged hadron and a single residual photon with a probability 
P 7 greater than 0.2 are rejected. 

The K* (decaying into K\ 7r) component is reduced by a cut on the HCAL energy 
deposit (-Ehcai) and the azimuthal offset (<50) between the track impact and the energy 
weighted barycentre of the energy deposit in HCAL 
following variables: 



25fl . The cut is applied on the 



-Ehcal — Ph 
Oh 

C 1 0barycentre 



(7) 



-'track impact | v"/ 

where Ph is the hadron momentum and the corresponding HCAL energy resolution. 
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Figure 22: The energy dependence of the ratio of data to Monte Carlo efficiencies for the 
additional cuts in the pion channel in the 1991-1995 data. 

The value of the parameter £ is +1 if 0barycentre is in the direction of the charged hadron 
bending, and —1 otherwise. The following region in the (Se ,84>) plane is cut : 

5 E > 
5<p < -1° 
5<p < (ffe-3)° 

In order to verify that those supplementary cuts do not introduce biases into the 
measurement of the r polarisation, the ratio of the efficiencies for each of the three 
cuts on data and on Monte Carlo, edata/^MC, is studied as a function of the pion 
momentum (Fig. [22|). The computation of these efficiencies and of their statistical 
errors takes into account the non-pion background estimated by Monte Carlo. The 
small bias introduced by these cuts is evaluated by a linear fit and corrected in the final 
polarisation measurement. The corresponding statistical error is taken as a systematic 
uncertainty. 

For the p channel, a cut on the invariant mass is required, m p > 0.4 GeV/c 2 , to 
reduce the feedthrough from the h mode. For the a\ channel, the invariant mass is 
required to be less than the r mass. In the case of a decay into three charged hadrons 
(3h), three good tracks are required in the hemisphere. 

The channel identification efficiencies, the r background, and the non-r background 
contributions as measured in the data (except for qq) are shown in Table |7|. 
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5.3 The r direction 



5.3.1 Kinematic reconstruction of the r direction 

For events where both r's decay into hadrons (called hadronic events), it is possible to 
reconstruct the t + t~ direction, thereby enhancing the sensitivity of the polarisation 
measurement in the p and a\ channels as seen in Section |2|. This method is used here 
for the first time. 

In hadronic events, although the momentum of each neutrino is unknown, the r 
direction can be determined up to a twofold ambiguity. The two solutions lie at the 
intersection of the cones having for axis the hadronic direction in each hemisphere and 
an opening angle <fio computed from the measured hadron momentum Ph and energy 
Eh, assuming the pion mass for each individual hadron, 

2E T E h — m T 2 — m h 2 
cos 0o = 2P~P h ' ^ 

where P T and E T are the energy and the momentum of the parent r assumed to be 
produced at the beam energy, and nth the hadronic mass. 

A procedure is used for events where the cones do not intersect, mainly because 
of detector resolution effects. The hadron momenta are allowed to fluctuate within 
the expected resolutions of charged and neutral pions. The r direction is taken as the 
average of the found solutions among the random trials, provided at least 10 solutions 
are found for 500 trials. This procedure attributes a r direction to half of the events 
where the two cones are not initially intersecting. 

As a result, 80% of the hadronic events are available for the polarisation analysis 
using the r direction. The remaining 20% without a reconstructed direction are 
analysed in the standard way. On the whole 52% of the r —>■ pv decays benefit from 
this improvement. 



5.3.2 Handling the two r directions 

Since the typical r decay length is 2 mm at LEP I, the precise determination of 
the secondary charged tracks brings some information allowing in principle to lift the 
twofold ambiguity and to choose the actual r direction. This procedure stems from 
the 3-dimensional method developed for the measurement of the r lifetime 0. In 
practice one can only separate the two solutions on a statistical basis with the help of 
an estimator, Dh, and assign to each of the two directions a probability P\^ to be the 
true one. This procedure is explained in Appendix A, while more details can be found 



in Ref. 26 . 



Reference distributions for the estimator are set up with Monte Carlo for various 
event configurations (vr-7r, p-p, p-n, etc.) and probabilities P\ 2 for each of the two 
found t directions to be the true one are obtained. 

If the highest probability is retained, it is found by simulation that the closest of the 
two directions to the true one is chosen in 65% of the cases for all channels combined. 
A more relevant way to characterize the procedure is to consider the precision achieved 
on the r direction. While the mean value of the angle between the two directions is 
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24 mrad, the angle between the direction with the higher probability and the true one 
is 12 mrad on average. The corresponding value for the direction closest to the true 
one is 7 mrad (dominated by resolution effects) and a value of 15 mrad is obtained if 
a random choice between the two solutions is made. In practice, the two solutions are 
retained, with their respective probabilities. 



5.3.3 Analysing the polarisation with the r direction 

The polarisation is analysed in each hemisphere using the proper optimal variables, 
ui t 2, calculated for the observed decay for each choice of direction using the description 
of the hadronic system given in Ref . || . Both directions are entered into the expected 
decay distribution depending on P T : 

W = P 1 F(u 1 )(l + P T oo 1 ) + P 2 F(lu 2 )(1 + P t lu 2 ) 

= F{l + P T u) (10) 

with the new optimal observable 

= PiFQjiVi + P 2 F(uj 2 )uj 2 
U P 1 F(lu 1 ) + P 2 F(u 2 ) ' 1 j 

The ideal sensitivity given in Table [T] is naturally degraded by detector resolution 
effects and mostly by the imperfect determination of the r direction. The expected 
sensitivities are given in Table |8| from a Monte Carlo study. The gain in sensitivity 
using the calculated probabilities instead of affecting each direction with an a priori 
50% weight is quite small, but it provides the ultimate gain in precision achievable 
within the ALEPH detector capabilities. 

For events where one of the r's decays leptonically and for which the information on 
the r direction cannot be retrieved, the acollinearity angle between the decay products 
in the two hemispheres contains additional information as compared to the separate 
analysis of single decays. This method B |23j is also used in this analysis. 



Table 8: Sensitivities obtained in the p and the a\ channels, for various methods using 
the true or the reconstructed direction. Only events with two reconstructed directions are 
considered here. 



method 


p-had 


ai(37r)-had 


ai(n2ir )-had 


co> truc , without f 
u truc , with t 


0.48 ±0.01 
0.57 ±0.01 


0.41 ±0.01 
0.57 ±0.01 


0.41 ±0.01 
0.56 ±0.01 


cu vec , without t 
uj rec , with Pi = P 2 = 0.5 
u; rec , with P 1>2 


0.46 ±0.01 
0.51 ±0.01 
0.51 ±0.01 


0.40 ±0.01 
0.47 ±0.01 
0.48 ±0.01 


0.37 ±0.01 
0.42 ±0.01 
0.43 ±0.01 



5.4 Summary of systematic effects 

Possible sources of systematic effects come from the r selection, particle identification, 
and 7T° reconstruction. The dominant systematic uncertainties are due to photon 
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Table 9: Summary of systematic uncertainties (%) for A T and A e with the r direction 
method. 



Systematic effect 


h 


1 

P 


3h 


h2ir° 


e 


n 


acol 


eff. h — *■ h id. 


0.17 


0.06 


- 


0.06 


0.20 


0.35 


0.01 


misid. (e, fi) — > /i 


0.24 


0.05 


- 


0.09 


0.13 


0.25 


0.57 


rr selection 


0.13 


0.03 


0.01 


0.01 


0.03 


0.04 


- 


r BR and background 


0.04 


0.05 


0.03 


0.09 


0.01 


0.02 


0.02 


tracking 


0.08 


0.07 


0.22 


- 


- 


0.21 


0.30 


7-reconstruction 


- 


0.22 


0.29 


0.66 


- 


- 


- 


7T -reconstruction 


0.11 


0.29 


0.68 


0.62 








fake photons 


0.31 


0.17 


0.28 


0.75 




- 


- 


ECAL scale 




0.20 


0.33 


0.63 


0.15 






ECAL + HCAL cut 


0.22 














modelling 






0.68 


0.68 








non-r background 


0.24 


0.16 


0.07 


0.05 


0.73 


0.50 


0.60 


r MC statistics 


0.34 


0.30 


0.61 


0.77 


0.73 


0.80 


1.44 


TOTAL 


0.66 


0.57 


1.30 


1.70 


1.07 


1.06 


1.69 


Systematic effect 


h 


A e 
P 


3h 


/i2tt° 


e 


P 


acol 


tracking 

non-r background 
modelling 


0.04 
0.13 


0.08 


0.02 
0.40 


0.07 
0.40 


1.23 


0.05 
0.24 


0.24 


TOTAL 


0.13 


0.08 


0.40 


0.41 


1.23 


0.24 


0.24 



identification and 7r° reconstruction, the non-r background, and, in the case of the a\ 
modes, the dynamics of the decay. Table ^ shows the summary of the main components 
of the systematic uncertainties for A T and A e . 

The particle identification efficiency matrix is measured as a function of momentum 
using data samples of kinematically identified particles, as described in Section f>. 2. 2 . 
The procedure yields both the relevant efficiencies and their associated uncertainties 
which are dominated by the statistics of the control samples. The effect of this 
systematic uncertainty on the polarisation measurement is significant only for the e 
or /i channels. 

To study the effect of the r pair event selection, the difference in the selection 
efficiency between data and the Monte Carlo is analysed for each cut, as discussed in 



detail in Refs. |[L7| , |19| . The systematic error on the selection has a significant effect 
only on the h channel where the nearly flat momentum distribution is distorted at low 
values by the 77 cuts and for large momenta by the Bhabha and //-pair cuts. 

The branching ratios used in the Monte Carlo are corrected using the measured 
values of ALEPH |I7 , |i~9f and the uncertainty in the amount of feedthrough in a given 
r decay channel from the other modes is obtained by varying the branching ratios 
within the errors given in Refs. [IT], taking into account the correlations between 
them. The impact on the r polarisation measurement is negligible. 
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The tracking resolution affects only the a,\ decay since three reconstructed charged 
hadrons are required. In this decay the tracks are generally close to each other and 
therefore partially overlap. A reconstruction problem could occur if the two same-sign 
tracks have nearly the same transverse momentum. Data and Monte Carlo distributions 
of the angle between the two tracks at the vertex are in fact in agreement within 
statistics and a systematic uncertainty for the polarisation is derived. The systematic 
error on the momentum calibration is also estimated with the largest effect observed 
on the fj, channel. 

A complete investigation of the observables entering photon and ir° reconstruction 
is undertaken. For each variable a comparison between data and Monte Carlo is 
performed in order to search for possible systematic deviations. These differences are 
then injected in the complete reconstruction procedure, affecting both the number 
of events in a given channel and the corresponding uj distributions. Each data- 
Monte Carlo comparison yields a systematic uncertainty limited by the statistics of 
the data. As discussed in full detail in Ref. fl7| , these studies include the effect of 
the clusterization algorithm, the photon probabilities, the threshold energy cut for low 
energy photons or the minimal distance between a charged track and a photon in the 
electromagnetic calorimeter. Similarly, the effect of the cuts on the ir° estimators is 
also determined. All the contributions to the total systematic uncertainties for this 
source are added in quadrature. 

Another source of uncertainty in the photon and tt° treatment is due to the 
problem of fake photons generated by hadron interactions in ECAL or electromagnetic 
shower fluctuations. This source of photon candidates is underestimated in the 
Monte Carlo simulation with respect to the data. This deficit produces a bias on 
the polarisation measurement and needs to be corrected for with the systematic 
uncertainty carefully estimated. In general, the effect has been studied by removing 
fake photons in the simulation and repeating the analysis: decay classification and 
polarisation determination in all channels. This method establishes the sensitivity of 
the measurement to the problem. Then the fake photon deficit in the simulation is 
measured in all the needed hadronic channels by fitting probability distributions in 
data and Monte Carlo as discussed in Section |5.2.4| and the corresponding systematic 
uncertainties are deduced. 

A more direct method is used in the h channel for the events with a residual 
photon where a tighter probability cut was applied to reject more p background (see 
Section |5.2.5|) . Since the p contribution can be well estimated by Monte Carlo, it is 



possible to use the results of this cut both on data and Monte Carlo to estimate the 
relative effect of fake photons as a function of momentum. Any momentum dependence 
of the Monte Carlo deficit would produce a bias for the polarisation. Such a bias is 
measured (Fig. |22|) and the statistical uncertainty in the linear momentum fit is taken 
as systematic error. This procedure in the h channel does not take into account the 
loss of events with multiple fake photons entering the definition of 7r° candidates. A 
special study was made to estimate this effect as the simulation tends to underestimate 
the fake photon multiplicity, with the result that the effect on the polarisation is small 
compared to that of residual photons. 

Since some additional calorimetric cuts are introduced in order to further reduce 
the level of the p and K* background in the hadron channel, it is important to study 
their effect on the hadron momentum distribution. The results have been presented in 



37 



Section |5.2.5| and the corresponding systematic uncertainties derived. 

Finally, the effect on the A T and A e measurements from the different sources of non- 
r background has been studied in detail directly on the data with the method presented 
in Section |5.2.3| . This procedure allows the energy distribution of the remaining Bhabha 
events in the selected sample to be investigated and a direct measurement of their 
cos# distribution to be made. The latter information is particularly relevant for the 
A e determination because of the asymmetric Bhabha contamination. Using as input 
the distributions derived from data with their statistical and systematic uncertainties, 
their impact on A e and A T is assessed. It is important to note that the dominant 
statistical uncertainty affects both the normalization of the background and the shape 
of the uj and cos 9 distributions. The most seriously affected channel is clearly that of 
the electron, both for A T and A e . 



5.5 Results 

Table |10| gives the results for A e and A T in the seven channels considered in this 
analysis. The statistical and systematic uncertainties are given separately. The last 
line presents the combination of these results taking into account their correlations, 
including the polarisation correlation between the two hemispheres in an event. 



Figures |23], |24], and ^5] show the u distribution for the p and a\ decay modes. The 
u distributions for the decays where it was possible to reconstruct the r direction are 
also shown separately. 



Table 10: A e and A T results with statistical and systematic uncertainties for the 1990-1995 
data with the r direction method. 



Channel 


Ar (%) 


A e (%) 


hadron 


15.49 ± 1.01 ±0.66 


17.36 ± 1.35 ±0.13 


rho 


13.71 ±0.79 ±0.57 


15.04 ± 1.06 ±0.08 


al(3h) 


15.01 ± 1.55 ± 1.30 


15.78 ±2.07 ±0.40 


al(h2vr°) 


15.94 ± 1.73 ± 1.70 


12.65 ±2.31 ±0.41 


electron 


14.58 ±2.18 ±1.07 


16.67 ± 2.92 ± 1.23 


muon 


14.45 ± 2.13 ± 1.06 


12.05 ±2.78 ±0.24 


acollinearity 


13.34 ± 3.83 ± 1.80 


19.41 ±5.02 ±0.24 


Combined 


14.58 ±0.53 ±0.39 


15.50 ±0.71 ±0.11 



6 Combined results 
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(a) Events without r direction. 
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Figure 23: Distributions of the u observable used in the polarisation fit of the p channel. 
The data are shown by points with statistical errors bars. The dotted and dash-dotted lines 
corresponds to the contributions of left- and right-handed r's respectively, as fitted in the 
data. The shaded area shows the non-r background contribution. The solid line corresponds 
to the sum of all simulated contributions. 
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(b) Events with r direction. 
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(c) Standard uj for all events. 



Figure 24: Distributions of the uj observable used in the polarisation fit of the a\ — > 3^^ 
channel. The data are shown by points with statistical errors bars. The dotted and dash- 
dotted lines corresponds to the contributions of left- and right-handed t's respectively, as 
fitted in the data. The shaded area shows the non-r background contribution. The solid line 
corresponds to the sum of all simulated contributions. 
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Figure 25: Distributions of the u observable used in the polarisation fit of the oi — > 7r ± 27r° 
channel. The data are shown by points with statistical errors bars. The dotted and dash- 
dotted lines corresponds to the contributions of left- and right-handed r's respectively, as 
fitted in the data. The shaded area shows the non-r background contribution. The solid line 
corresponds to the sum of all simulated contributions. 
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6.1 Standard model corrections 

The A T and A e parameters are obtained by fitting the integral, in each cos 8 bin, of 
the function (|l|) to the measured polarisation. To obtain the combinations of effective 
couplings (Aj) related to sin 2 6^y, corrections have to be applied which take care of 
the fact that the polarisation was not measured at the Z pole and that there are 
radiative corrections and Z-7 interference effects. These corrections are computed by 
the ZFITTER program [f|] in the following way. 

For a set of the standard model parameters (Z mass, top mass, Higgs mass, a s ), 
ZFITTER provides the value of the effective couplings in the framework of the standard 
model. Cuts are defined on acollinearity and momentum which reproduce the ones used 
in the analysis. Then the Afb asymmetry, and the polarisation in each cos 8 bin are 
calculated at the given energy. The same fitting procedure as for the data is applied 
and the difference between the effective couplings and the fitted parameters is taken 
as the correction. 

The study of the stability of the corrections against the variation of the input 
parameters within their errors shows that the related error is negligible compared to 
the systematic errors involved in the analyses. The same is true for the values chosen 
for the cuts. With this procedure, the inaccuracy of ([!]) in describing the polarisation 
does not introduce any bias but may simply reduce very slightly the sensitivity. The 
level of these corrections is 0.04% with an uncertainty of 0.01%. It has been checked 
that using as input the measured values of Afb gives the same results. 

6.2 Consistency and combination of the results 

To investigate the consistency of the two analyses, both statistical and systematic 
effects have to be taken into account. 

As a result of the different selection procedures, there is a significant number of 
events which belong to only one of the two samples. The fractions of unshared events 
in the different channels, averaged over the full angular distribution, are 24% (tt), 
27% (p), 24% (a x -> 3vr), 47% (a x -> tt2tt°), 24% (e), and 13% (//). Typically, these 
fractions are by a factor 1.6 larger in the small-angle region (| cos^l > 0.7) than in 
the central part. The distribution of the events has to be taken into account in the 
evaluation of the consistency and in the combination of the measurements since this 
small-angle region has a significant contribution to the A e determination. Another 
statistical effect comes from the channels where the r direction is taken into account 
since different information is used in the two analyses. Finally, the Monte Carlo samples 
do not completely overlap. 

The systematic effects of the two methods have in common the detector behaviour 
and the basic event reconstruction; however the analysis tools used are largely 
independent. Furthermore, the procedures developed to measure efficiencies, and 
various corrections are different in the two analyses and thus a part of the estimated 
systematic uncertainties is uncorrelated. 

The A T values given by the two analyses - for each channel and globally - are 
in agreement within the estimated uncertainties. The combination of all channels, 
gives AA T = (0.14 ± 0.52)%, where the statistical and systematic part of the error 
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are comparable. The agreement is not as good in the case of the A e measurement 
where the deviation between the two analyses, mainly due to the hadron channel, is 
AA e = (0.92 ± 0.48)%. However, the systematic uncertainties on A e , dominated by 
the effect of the non-r background (chiefly Bhabha events), are smaller by a factor 
of 5-10 than the statistical errors, and the many cross-checks performed show that 
the background levels and distributions are known within the quoted uncertainties. 
This lead to the conclusion that the observed difference is the result of a statistical 
fluctuation. 



6.3 Final results 

The consistency of the two sets of results having been checked, they are combined 
to obtain the final results. Since the two analyses have comparable statistical and 
systematic uncertainties, the combination procedure is not a critical one. It was chosen 
to average the results using equal weights. For the computation of the final errors the 
statistical correlation between the two samples is taken into account together with the 
systematic uncertainty correlations. 
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Figure 26: Polarisation dependence on cos 6 for the LEP I data. The curves corresponding to 
Eq. [l] for parameters given by fits with (dashed line) and without (dotted line) the universality 
constraint are superimposed. 



The results are given in Table [IT] for all the analysed modes and in Table [T2| for the 
global values. Figure p6| shows the average tau polarisation as a function of cos#. 
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7 Discussion 



The A T and A e parameters are measured by means of similar methods by the three 
other LEP collaborations [27-29]. The A e asymmetry is also measured by the SLD 
collaboration PU| using the beam polarisation. 
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Figure 27: (a) Comparison of the A e and A T measurements by the ALEPH (A), DELPHI 
(D), L3 (L), and OPAL (O) collaborations. The ellipses are one standard error contours 
(39% CL). The horizontal lines represent the SLD (S) A e measurement through the left-right 
asymmetry of the cross-sections plus and minus one standard deviation, (b) Comparison of 
ALEPH A T , A e , and A® b measurements. The ellipse is the one standard error contour. The 

hyperbolas and horizontal lines represent the A°j b and Aj b 
and minus one standard error). 



"" ~ 1 ' n,e measurements respectively (plus 



All the results are shown in Fig. 
be observed. 



(a) where a reasonable general agreement can 



The A T , A e , and A°j b asymmetries at the Z pole are related by the equations: 



A 



0,e 



-{A e 



A 



O.r 



:i2) 



The comparison of the present measurements of A e and A T with the forward-backward 



asymmetries measured by the ALEPH collaboration [11] is presented in Fig. |27| (b). 



Assuming universality, the leptonic forward-backward asymmetry measured by 

0,1 



0.0173 ±0.0016 |lll, to be 



ALEPH, using all the leptonic Z decay channels, is Aj b 

compared with J(A-r) 2 = 0.0163 ± 0.0010. 

The information given by the r polarisation measurement can be used, together 
with the forward-backward asymmetries and the partial widths of the Z decays into 



lepton pairs [pdj| , to determine the Z leptonic couplings. The results obtained by the 



procedure of Ref. [11] are given in Fig. and Table 
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Table 11: Combined A T and A e results with statistical and systematic uncertainties for the 
1990-1995 data. The pion inclusive channel is highly correlated with the hadronic channels. 



Channel A T (%) A e (%) 



hadron 


15.35 


± 


0.96 


± 


0.50 


16.32 


± 


1.27 


± 


0.09 


rho 


13.75 


± 


0.78 


± 


0.42 


14.85 


± 


1.03 


± 


0.06 


al(3h) 


14.89 


± 


1.50 


± 


1.05 


14.68 


± 


1.98 


± 


0.40 


al(h27T°) 


16.14 


± 


1.74 


± 


1.44 


14.13 


± 


2.30 


± 


0.42 


electron 


14.11 


± 


2.18 


± 


0.84 


15.38 


± 


2.92 


± 


0.77 


muon 


14.05 


± 


2.07 


± 


0.82 


11.91 


± 


2.71 


± 


0.17 


acollinearity 


13.34 


± 


3.83 


± 


1.80 


19.41 


± 


5.02 


± 


0.24 


pion inclusive 


14.93 


± 


0.83 


± 


0.87 


14.91 


± 


1.11 


± 


0.17 



Table 12: Combined A e , A T and A e - T results with statistical and systematic uncertainties 
for the 1990-1995 data. 



A T (%) 14.51 ±0.52 ±0.29 

A e (%) 15.04 ±0.68 ±0.08 
Ae- T (%) 14.74 ±0.41 ±0.18 



Table 13: Effective lepton vector and axial vector couplings. 



9v 


-0.0384 ±0.0015 


9a 


-0.50171 ±0.00056 


fv 


-0.0377 ±0.0055 


9\ 


-0.50043 ± 0.00094 


9v 


-0.0367 ±0.0015 


9 t a 


-0.50219 ±0.00093 


9 l v 


-0.03747 ±0.00096 


9 1 a 


-0.50155 ±0.00045 
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Figure 28: Effective lepton couplings. The ellipses are one standard deviation contours (39% 
CL). The shaded area indicates the Standard Model expectation for M t = 174 ± 5 Gev/c 2 
and 90 < M H (GeV/c 2 ) < 1000; the vertical arrow shows the change if the electromagnetic 
coupling constant is varied within its error. 

Comparing these numbers with the values of the couplings presented in Ref. |TT 



shows the important weight of the r polarisation measurement in the Z leptonic 
coupling determination. 



8 Conclusion 

From the A T and A e values obtained through the r polarisation measurement: 

Ar = 0.1451 ± 0.0059 , A e = 0.1504 ± 0.0068 , 
the following ratios of the effective couplings are computed: 



9v/g T A 

9v/9 e A 



The universal value 



0.0729 ± 0.0030 
0.0756 ± 0.0035 

0.964 ±0.060 . 



A e - T = 0.1474 ±0.0045 



9v/9a 
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translates into a determination of the effective weak mixing angle 



sin 2 9f r = 0.231 47 ± 0.00057 . 



(13) 



Using both the r polarisation and the leptonic forward-backward asymmetries 
measured by ALEPH, the value 



is obtained. 
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sin 2 ef r = 0.23130 ± 0.00048 



(14) 
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A Tau direction estimator 




Figure 29: Decay of a r pair. 



A.l The ideal case 

First consider the case of a r event where both r's decay into a single pion with 
momentum vector measured perfectly (Fig. |29[). Let To be the interaction point and 
H± points on each reconstructed track, in the vicinity of the decay points. The vectors 
/i ± and f ± are unit vectors along the 7r ± and r ± momenta respectively, l± are the r ± 
decay lengths and 5± the distances of H ± to the r ± decay points. From f + = — r~ 
and the relation (Fig. |2S|) 

of + / ± f ± = on± + sji* , 

where O is an arbitrary origin, the equation 

It = SJr - 5+h + + H^lL (15) 

follows, where / = /+ + /_ and f = f ~. 

Taking for H ± the end points of the vector D h defined ||[] as the minimum approach 
vector between the two reconstructed tracks, the previous equation becomes 

lr = lA h + D h , (16) 

where the vector is orthogonal to the hadronic plane spanned by h + and h~ and 
the vector Ah lies in this plane. 

The angles ± (Fig. |2S]) are given by Eq. [5], 

, 2E T ±E h ± - ml- m 2 h± 
cos^ = T - 2- , 

1 A similar analysis of the ideal case is given in Ref. [pl| . 
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hadronic plane 



ambiguity plane 



Figure 30: Definition of vectors Ah, Dh, and N. 

consequence of the r mass constraint, where E, p, and m are energy, momentum, and 
mass of r's and hadrons (vr's). Accordingly, the two kinematically found r~ directions, 
defined by the unit vectors fx and f 2 , are the intersections of two cones having — h + 
and h~ for axes and ^ for opening angles. The vector f\ + f 2 lies then in the hadronic 
plane and is, therefore, orthogonal to Dh (Fig. |3(]). This implies 

D h -n = - D h -f 2 . (17) 

The vector Ah, projection of f on the hadronic plane, is the internal bisector of T\ and 
f 2 . These properties, with Eq. [16|, lead to 



\D h 



1 - lAfcl 



since I is, by definition, a positive length. The positivity of I and Equations [L6| and ^7 
imply that the physical r direction is determined by the relation 

D h -t truc >0. (18) 

For hemispheres with a 3-prong decay, the vector h ± is defined as equal to the sum 
of the three momenta and originating from the 3-track common vertex. 



A. 2 The general case 

When 7r°'s are produced, the vector Ah, which is kinematically determined, is still 
measurable, but the exact reconstruction of Dh given in the previous section is not 
possible and an approximation must be made. 
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The vectors are now unit vectors along the sums of the hadron momenta in 
each hemisphere. Two vectors, which make with Ah an orthogonal system, can be 
unambiguously measured: 

h± = h~ x h + , 
which is normal to the hadronic plane, and 

N = h± x t = h + cos (f)~ + h~ cos <p + , 

which is normal to the ambiguity plane spanned by f\ and fi- In terms of h± and N, 
the analytical expressions of Ah and Dh are 

A h = Nxh ± /h 2 ± , 

D h = h ± (M-h ± )/h 2 ± . 

The vector Aft is equal to H+H_, where H± are arbitrary points on the lines 
drawn along the vectors h from the r ± decay points, lines which coincide with the 
reconstructed tracks in the ideal case of the previous section. 

The charged tracks reconstructed in each hemisphere are used to define an 
approximation of Dh- Let c be unit vectors along the momenta of the charged 
particles, C± points on the char ged tracks in the neighbourhood of the decay points, 
C± = C~ X c , and A& = C^(L. If the vectors c± and h± are collinear, the space 
planes parallel to the hadronic plane and containing the r ± decay points contain also 
the charged tracks. The vector Dh is then equal to the minimum approach vector 
between the charged tracks. In the most general case, an equation similar to Eq. 
can be written |J: 

If = IA C + D c , (19) 

where the vector D c , whose length is proportional to the minimum distance between 
the charged tracks, is collinear to h±: 

A c = Nxc ± /(h±- c ± ) , 

D c = h ± (K6- c ± ) I (h± ■ c ± ) . 

As the relation (|T6|) is still true, the expression of the Dh vector can be written as : 

D h = D c + l (A c - A h ) . (20) 

In the equation (pQ[) , the vectors D c , A c , and Ah are measurable, but the decay 
length remains unknown. Taking the mean r decay length sum I = 4.6 mm as a value 
for I, an approximation of Dh, denoted Dh aH , is obtained. 

The vector A c lies in the ambiguity plane and, since both Dh and D c are collinear to 
h±, the projection of IA C on Ah is lAh, analytically A c - Ah = N 2 /h 2 ± = A\. Therefore, 
the vector D hcB is proportional to D h , irrespectively of the value of I. As in the ideal 

case, 

D Ks ■ fi = ~Dh cB ■ f 2 , (21) 

and the r direction giving the positive value for Dh eS ■ f can be selected. 

Because of the approximation made when 7r°'s are present and resolution effects, 
this selected direction is only correct on a statistical basis and probabilities, computed 
from reference distributions of Dh cS ■ f, are assigned to the two solutions. 
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